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Abstract. We study a new Radon-like transform that averages 
projected p-forms in R" over affine (n — fc)-spaces. We then prove 
an expHcit inversion formula for our transform on the space of 
rapidly-decaying smooth p-forms. Our transform differs from the 
one in [GGS]. Moreover, if it can be extended to a somewhat larger 
space of p- forms, our inversion formula will allow the synthesis of 
any rapidly-decaying smooth p-form on R" as a (continuous) su- 
perposition of pullbacks from p-forms on fc-dimensional subspaces. 
In turn, such synthesis implies an explicit formula (which we de- 
rive) for reconstructing compactly supported currents in R" (e.g., 
compact oriented fc-dimensional subvarieties) from their oriented 
projections onto fc-planes. 



1. INTRODUCTION 

Question: Do the projections of an oriented p- dimensional object in 
R" into /c- dimensional subspaces — cancellation allowed — determine it? 

Question: Can one use its projections to explicitly reconstruct the 
object? 

Here, these questions prompt us to define a new integral transform 
of "Radon" type for differential p-forms, and to derive a formula that 
inverts it on the space of rapidly-decaying p-forms. Our inversion for- 
mula doesn't quite settle the questions above, but it comes very close: 
We believe it shows exactly what their answers have to be. 

By no means are we first to construct a tomography of p-forms: In 
the late 1960's, Gelfand, Graev and Shapiro did so in [GGS], and their 
work was followed by a series of related papers (e.g., [GGG]). But 
our theory complements theirs, and unlike it, bears directly on the 
questions above. 
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To understand our transform, and see how it differs from the one in 
[GGS] , recall that any a linear /c-plane P , together with its orthogonal 
complement P' , splits R" as the direct sum R" = P ® P' . The 
exterior algebras f\^P and then both form subspaces of the 

larger exterior algebra ■ Note, however, that their sum does not 

generally span A*!^" ■ 

Now consider a continuous p-form a in R" , and suppose it decays 
fast enough to make it integrable over any affine (n — A;)-plane. Our 
transform assigns, to each ^ G P, a value TZko{P.^) G /\pP , by 
projecting a{x) into /\pP at each x E P' , and then integrating along 
P' . We emphasize that TZka{P,S) thus belongs to f\pP ■ 

Like ours, the transform in [GGS] integrates along P' . But it does 
so after projecting a{x) into f\^P' — not /\^P — at each x E P' . 
Accordingly, their transform assigns, to each ^ e P , a value in /\pP' 
instead of Ap-^ ■ Though procedurally similar to ours, this process 
encodes very different information about the input p-form a . 

For instance, when a is a 1-form in R^ , our 2-plane transform (which 
integrates along lines perpendicular to each linear 2-plane) is invert- 
ible, in the sense that we can recover any rapidly decaying 1-form a 
from its transform 7^2« • The transform in [GGS], contrastingly, is not 
invertible in that case. For it annihilates all exact 1-forms. Indeed, its 
integration step computes standard line integrals, and the line integral 
of a rapidly decaying exact form must vanish. 

Further, the low dimensions of the example just given oversimplify 
the distinction between our transform and that of [GGS]. For when 
k = n — 1 = 2 , and p = 1 , the inclusion 

becomes an isomorphism. In that case, adding our transform to the 
one in [GGS] gives the classical (trivially vectorized) Radon transform. 
Roughly speaking, that puts any one of these transforms in the span of 
the other two, in which case our work produces nothing fundamentally 
new. But in general — in fact, whenever k > p > 1 — the inclusion 
above is proper, and this makes our transform truly inequivalent to 
that of [GGS]. 
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To explain why our transform and inversion formula are not just new, 
but also interesting, we revisit our opening questions. The oriented 
p-dimensional objects we have in mind there are compactly supported 
p- dimensional currents, in the sense of deRham. As we explain in §6, 
they belong to the dual of the space of smooth differential p-forms, and 
they project onto /c-dimensional subspaccs in a simple, natural way. 
Their relationship to p-forms makes the problem of reconstructing a 
current from its projections equivalent to that of writing any rapidly 
decaying p-form as a superposition of forms that are k-planar, that 
is, puUbacks of p-forms from /c-dimensional subspaces via orthogonal 
projection. Formally, one gets precisely the right type of superposition 
formula from our inversion result, and when p = (where our whole 
theory reduces to that of the classical scalar Radon transform) this 
formula is known to hold. We feel sure that it holds for p > too, 
and state the expected facts, along with the answers they give to our 
opening questions, as explicit conjectures after proving our inversion 
result (Theorem 6.1). We hope that our work tempts someone having 
more skill with the required analytical tools to tie up these loose ends. 

Meanwhile, discovering the appropriate tools, deriving the correct in- 
version formula, and proving its validity on the smaller space of rapidly 
decaying p-forms will occupy us here. Before beginning, we briefly pre- 
view our plan of attack. 

Section 2 sets preliminary concepts and notation dealing mainly with 
exterior algebra and the Fourier transform. The material is routine, ex- 
cept possibly for Propositions 2.8 and 2.9, which show how the Fourier 
transform interwines d , S , dS , and Sd with certain operations in the 
exterior algebra. These elegant formulae are easy, but we have not seen 
them elsewhere, and they combine to intertwine the half-laplacians d6 
and 6d with the projections we call 11 and $ in . We later ex- 

ploit these facts to derive Corollary 6.2, a simpler "differential" state- 
ment of our inversion formula valid for even codimension. 

In §3, we define our transform TZk and its dual TZl . We work out a 
couple examples along the way to clarify definitions and suggest the 
typical behavior of these operators. 
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In §4, we start the core work of our paper by deriving an explicit 
formula for the composition TZl^oTZk as a convolution (Theorem 4.4). 
This type of result — and its utility for inverting a Radon transform — 
were first noted by B. Fuglede in [Pu]. Indeed, our formula reduces to 
Fuglede's result when the p-forms on which it operates are mere scalar 
functions (p = 0). But the precise generalization for p > requires 
a more sophisticated analysis focused on the action of the orthogonal 
group 0(n) on the exterior algebra A*!^"- The resulting convolution 
formula thus forms the first significant contribution of our paper. 

The other really new element is our computation of the precise Fourier 
transform of TZloTZ^, obtained as Theorem 5.1 in §5. Roughly speak- 
ing, we carry it out by expanding the convolution kernel from §4 as 
a sum of quotients of spherical harmonics by powers of the radial dis- 
tance r . We then use a well-known lemma from Stein [St] to compute 
the Fourier transforms of these quotients. 

We conclude in §6 by combining the convolution formula and our 
Fourier transform computation to invert TZk (Theorem 6.1). When 
the codimension n — k is even, our inversion formula simplifies: The 
intertwining rules at the end of §2 replace the Fourier transforms and 
projections there to simple linear combinations of the "half-laplacians" 
d6 and 6d . Corollary 6.2 states the precise result. 

Finally, our Inversion Theorem leads quickly to the formal superposi- 
tion result mentioned earlier: We state it as Conjecture 6.5, and then 
show that it would answer our opening questions. We present those 
answers too, as Conjecture 6.8. 

We completed this project with sabbatical support from Indiana Uni- 
versity, a Lady Davis fellowship at the Technion in Haifa, Israel, and 
extended hospitality from Stanford University. We deeply appreciate 
the assistance of all three institutions. 



2. Preliminaries 



For later reference, we record basic notation and background here, 
mostly about the exterior algebra, the Fourier transform, and their 
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interaction. The reader may prefer skipping to §3, and return here for 
clarification when the need arises. 

2.1. Subspaces and complements. Let Gn,k denote the compact 
Grassmann manifold of all A;-planes in R" . We always denote the 
orthogonal complement of a linear A;-plane P € Gn,k by P' , and we 
write k' for dim P' = codim P = n — k . 

2.2. Exterior algebra. We write A*^ ^^r the exterior algebra of 

a finite-dimensional vector space V . Chapter 1 of Federer's book [Fe] 
gives a careful (albeit terse) development of the exterior algebra from 
scratch. Here we review just a few relevant points. 

First, /\^V has the direct sum decomposition 



where /\fc^ denotes the subspace generated by fc-fold wedge products 
of the form f i A ^2 A • • • A , with Vi & V . We call the elements of 
/\f^V A;-vectors. 

Given an orthonormal basis ei,e2,...,e„ for V , we get a basis for 



where A runs over all possible increasing, length-A; subsequences of 
{1, 2, • • • ,n} . Declaring this basis orthonormal extends the inner prod- 
uct on V to an inner product (•, •) on /\^V 

We get a linear transformation on /\^V by left- wedging with any fixed 
vector V . We indicate the adjoint of this mapping (relative to the 
inner product above on /\^V) by "Vv" and define it by requiring 



Any linear map / : — > W between vector spaces induces a multilin- 
ear extension between their respective exterior algebras via the simple 




(2.2.2) 

for all a,l3e /\^V . 



{aVv, P) := {a, v A P) 



formula 



A ^2 A • • • A Vfe 



f{v^)Af{v2)A■■■Af{vk) • 
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Two examples particularly relevant for us arise with respect to a sub- 
space P C R" : 

(1) Inclusion: The inclusion f : P ^ R"- extends to an inclusion 
/\^P C A*R" ■ This lets us regard elements of /\^P as belong- 
ing to A*R" ■ We shall do so below routinely without further 
comment. 

(2) Projection: The orthogonal projection R" — > P is adjoint to 
the inclusion P C R" with respect to the usual dot product. 
It extends to a projection A*R" ~^ f\*P ^^^^ orthogonal 
for the inner product (•, •) defined above on A*R" • We often 
suffix a vertical bar to indicate this projection. For instance, 
we write the projection of a vector v on R"^ to P as v\p. 

Let Hx denote the hyperplane perpendicular to a non-zero x e R" . 
The sphtting R" — H^® H'^ induces a spht in the exterior algebra: 



We designate the projections onto the first and second factors here 
respectively by 



(2.2.4) n, : A.R'^ ^ A*^. and : A.R'^ ^ ^ A A.^. ■ 



The first of these is induced by the orthogonal projection R" — >■ H^. , 
as in item (2) above. We give it a special symbol nonetheless — it plays 
a key role later on. 

Finally, note that acts as the identity — and vanishes — on the 
degree-0 summand Ao-^" ' because the latter hes entirely in f\^Hr^ . 

2.3. Differential forms. For our purposes, a differential p-form on 

a subspace P C R" is a mapping (f) : P ^ ApR"* > where R"* denotes 
the dual of R'* . The dot product identifies R'* with R"* via 

{P1,P2, ■ ■ ■ ,Pn) ^ PidXi+P2dX2^ \-pndXn, 

and because of this isomorphism, we will frequently ignore the distinc- 
tion between vectors and covectors. We often regard maps from R" 
to ApR" as differential forms, for instance. 
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Call a differential form on R" A;-plancir if we can write it as the 
"puUback" of a differential form native to some /c-dimensional sub- 
space P C R" via orthogonal projection R"^ — > P . Using the notation 
introduced above, this simply means that for some P £ Gn,k a-nd all 
X e R" , 4> satisfies the following two conditions: 



For instance, the forms commonly denoted by dx + dy and xdx Ady 
are both 2-planar on R^ with P defined by z = . 

Contrastingly, neither zdxAdy nor dx A dy A dz are 2-planar. The 
former satisfies condition (2) for the plane z = 0, but not (1). The 
latter satisfies (1) but fails condition (2) for all P e G2,3 . 

Remark 2.4. The point of Conjecture 6.5 is that we may construct 
many g-forms in R" by averaging /c-planar g-forms over Gn,k ■ 

Next, recall that we take the exterior derivative of a differential form 
on R" using the formula 



Somewhat less familiar is the dcgrcc-dccrcasing divergence operator 
6 . Our notation (cf. 2.2.2) lets us to express it thus: 



Note that we omit the minus sign often used to define S so that 



making A denote the positive sum of pure second partials. In §2.7 we 
will see the Fourier transform relate S d and d5 to the projections 11 
and $ of (2.2.4) above. 

2.5. Mapping spaces. We will encounter tensor fields on R" , which, 
like differential forms, take values in a normed vectorspace V . To 



(1) ct>{x) = ct>{x\p) 

(2) <t>{x)^<t>{x)\^ 



(2.4.1) 




(2.4.2) 




(2.4.3) 
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measure decay rates of these mappings near infinity, we define the 
seminorm for any map F : — > V , and any s > , via 

(2.5.1) ||F|L . := hm sup \x\'\F{x)\ . 

Assigning a finite value to H-P'll^^) is shghtly more precise than saying 
"F is 0{\x\^^) ," which merely means < oo . Note that when 

F is locally integrable, we ensure summability on any proper subspace 
of — though not on R" itself — by assuming ||-F1|(„) < oo . 

We will say that F : R" — > V decays rapidly if H-P^Hj^^) < oo for 
every s > , and we define the space «S(R", V) of Schwartz TZ-fields 
as the set of all smooth F : — > V with rapidly decaying derivatives 
of all orders. To simplify notation, we will abbreviate 

S(W\ V, W) := 5(R", Rom{V, W)) , 

and when P C R" is a subspace, 

Sp{P) ■.= S(^P, ApP) (Schwartz p-forms) 
S{P) := So{P) = S{P, C) (Schwartz functions) 

2.6. Tempered {V, H^)-distributions. When V and W are inner- 
product spaces, we write (S*(R",V, ly) for the space of continuous 
linear functionals 5(R", V) ^ W , and call its constituent functionals 
tempered (V, W) distributions or Hom(V, iy)-valued distribu- 
tions, or, when context allows, simply tempered distributions on 
R'* . When applying such a distribution r to a Schwartz F-field F , 
we write t[F] , with square brackets to emphasize the distributional 
context. 

Certain distributions arise as locally integrable operator fields. An op- 
erator field is a mapping T on R" that assigns a linear transforma- 
tion Tx G Hom(V, W^) to each x G R" . When T satisfies IjT'llj^p-) < oo 
for some p > , and is also locally integrable, it represents a tempered 
Hom(y, iy)-valued distribution (which we still call T ) via integration: 

(2.6.1) T[F] := [ T{x) . F{x) dx , for all F e 5(R", V) . 
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Here and henceforth the "dot" operator in '^A.X" instructs us to 
feed the vector X e V to the operator A . This will save us from 
over-nesting parentheses. 

Equation (2.6.1) includes S{W,V,W) into S*(W,V,W). Since the 
inner product identifies V with its dual V* , it also identifies <S(R", V) 
with 5(R'^, V*) , so by (2.6.1) also includes <S(R", V) into 5*(R", V*) . 

As in the scalar setting, one can extend, to all tempered (V, l^)-distri- 

butions, many operators initially defined only on some subset that 
includes 5(R", V, W) . The classic example is differentiation: When 
T e »S(R", V, W) , integration by parts gives 



The right side above makes sense for any tempered {V, W^) -distribution 
T, assuming F G S(TV^,V), so the identity above serves to define 
dT/dxi for any T G 5*(R", V, W) . 

2.7. Fourier transform. Using the normalizing conventions of Stein 
[St], we define the Fourier transform of any integrable V^-field F on 
R** as the l^-field given by 



Any basic treatment (e.g. [Str]) will establish fundamental properties 
of the scalar [V = C) Fourier transform which extend trivially to our 
vector-valued setting. These include: 



(1) J" maps S(W,V) into itself. 

(2) is an isometry, and modulo a sign, inverts itself on 
S{W, V) . Specifically, 



The superscript "— " here signals composition with reflection 
through the origin, as in (7). 



(3) Given T e <S(R", V, W) and F e <S(R", V) , intertwines the 
convolution product T -k F , given by 




(2.7.1) 




F{x) e 



•^■^ dx . 



TT{F) = F 
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with functional operation, in the sense that 

T{Ti<F)= T{T).T{F) 

The definition of convolution doesn't actually require rapid de- 
cay of T and F , but when both are Schwartz, so is T -k F , 
which allows one to extend the Fourier convolution rule to a 
distributional setting. We prove general statements of this type 
in our Appendix (Observations 7.1 and 7.2). 

(4) is "self-adjoint" in the sense that when T e S{W, V, W) 
and F G 5(R"', V") , we have 

[ T{T) .F dx^ [ T. T{F) dx. 
JR" Jr" 

Note that by using distributional notation, we can express this 

last identity as 

(2.7.2) T{T)[F] = T[T{F)] , 

which extends the Fourier transform to all of <S*(R", V, W) . 

(5) Finally, recall the useful way that the T intertwines differenti- 
ation and monomial multiplication: 

For any F e »S(R", V) , we have the identities 

(BF \ B 
—J = 27ri HF) , ^^(F) = -27ri T {x, F) . 

These well-known identities follow directly from the definition of JF , 
but they imply less familiar rules for intertwining with the exterior 
derivative and divergence operators of (2.4.1) and (2.4.2): 

Proposition 2.8. For any (j) e 5'p(R") , we have 

T{d(j))= ■n:d\i\'^ :F{S(f))^ i7rJ^{(f))Vd\^\^ 

d:r(0) = -i7rJ^((i|a;p A0) (5 J^(0) = -iTrJT (0 v . 

Proof. Using d\x\'^ = '^^jXjdxj and similarly for the stated 

identities follow easily from (2.7.3) above, the adjoint identity (2.2.2) 
relating V to A , and the formulae for d and 5 given by (2.4.1) and 
(2.4.2). □ 
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Since d\x\'^ is dual to 2x , one easily combines the identities above 
with (2.2.2) to relate the "half-laplacians" d5 and 5d (2.4.3) to the 
projections 11 and $ as promised at the end of §2.3 above: 

Proposition 2.9. Let r :— \x\ denote the radial distance function on 
R" . Then we have the following identities on SpiYil^) : 

d5 = -Att'^ or'^^oT , 5 d = -47r2 JF-i o n o JP- 

and since d5 -\- 5 d — Is. , 

A = -47r^ JT-i o o JF . 

We don't know any reference for the first two identities here, even 
though the last is common knowledge. 



3. The transforms TZk and TZl. 

3.1. The transform. The transforms we describe in this section ex- 
change differential forms on R" with tensor fields on the canonical 
bundle Vn,k over Gri,k ■ We define the latter bundle via 

rn,fc:={(P,OeG„,fcxR" : eeP} . 

Definition 3.2 (The transform TZk)- Given a continuous differential 
p-form a on R** with ||q;||(^) < oo , we denote its transform by 

and define it as an integral: 



(3.2.1) 7^fca(P,0:= / + v) 

Jp' 



dr) . 

p 



Remark 3.3. When p = , so that a is just a scalar-valued function, 
we have a\p = a. In this case, TZk reduces to the classical A;'-plane 
transform (cf. [H, p. 28]). 

Remark 3.4. The definition above characterizes TZkCx as a mapping 
from r„ to A*-^" ' niaking it a section of the trivial A*-f^" bundle 
over Tn,k ■ Actually, TZua takes values in the sub-bundle whose fiber 
over (P, ^) is the puUback subspace /\^P C Ap R-"' ■ 
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Remcirk 3.5. The integral above does not signify oriented integration 
of forms in the Stokes-theoretic sense. For one thing, the integrand 
belongs pointwise to /\^P , hence vanishes on any /c'-tuple of vectors 
tangent to the domain of integration P' . Rather, the integral performs 
elementary vector integration of a /\pP-valued function against the 
Lebesgue measure on P' . The following example illustrates. 

Example 3.6. Suppose a e Ap^^" > A > , and define the p-form 
:= e~'^l^l^a . We compute TZkcf) using the well-known identity 

(3.6.1) / e-'\^\' dx ^ i^y'" , 

which implies, for any (P, ^) G r„ ^ , that 



I p' 
whence 



The input p-form above decays rapidly, and on each fiber of r„^fc , its 
transform TZk4> has the same property. The proposition below shows 
that this is no accident. 

Let |S"| denote the n-dimensional volume of the unit sphere in R"+^ . 

Proposition 3.7. Suppose a is p-form on R", and s>n. Then for 
any P e G„_fe , we have 

\\T^ka{P,-)\\^,_,,^<4\S''-'\ 

Proof. Without loss of generality, assume ||q;||(5) < oo for some s >n . 
Then for sufficiently large ^ and any 77 e P' , we have 



Ke + ^)I<2|HI(,) {i + \C\' + \v\'r^" 



hence also 



< 2H|(,)(l + |er)-^/^ / (^ + tK^) ^'^ drj 



1 + 
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Switching to polar coordinates and changing the radial variable r = I77I 
to p — rj a/TT"!^ , we now get 

Jo 



But 



POO POO 

J {i + p')-'^'p''-'dp < i + l {i + pT'^'p^'-'dp 



/OO 
p-'-^'-'dp , 

and —s + k' — l< —n + k' — 1 = —k — 1 < — 1 , because s > n . The 
last integral above is therefore bounded by < 1 . □ 

3.8. The dual transform TZ* . As in the classical Radon theory, our 
transform TZj. has a formal dual. In our situation, the dual sends 
continuous /\pR^ -valued fields on r„ ^ back to p-forms on R"^ . 

Recall that the transitive action of 0(n) determines a unique Haar 
probability measure on Gn,k ■ We compute all integrals over Gn,k 
using that measure. 

Definition 3.9. The dual transform TZl takes a continuous mapping 
P : rn,k — ^ Apf^" to p-form on R" given by 



dP , 

p 



Here Px denotes the orthogonal projection of x onto P . 
Note that we always have 



(3.9.1) niP = niP when (^(P, .) - /5(P, .)) 



= . 

p 



For each x G R" , the map P ^— > (P, Px) embeds Gn,k into r„ y^. . The 
compactness of Gn,k and continuity of /3 therefore guarantee existence 
of the integral defining TZl . 

Remark 3.10. Note too that for each fixed P G Gn.k ■ the integrand 
/3(P, Px)|p specifies a /c-p/anar p-form on R" . It follows that anything 
in the image of TZl is a superposition, over all P G Gn^k , of A;-planar 
p-forms. The statement of Conjecture 6.5 depends on this fact. 
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Example 3.11. Fixing an arbitary i> e , let us compute the dual 
transform TZ2P of the simple field P : r3^2 — Ai-^^^ given by 

By definition, 

r 1 /■ , , ,o 

du . 



nif3{x) = [ f3{P, Px) dP^^ [ e-l'^^^l'v 

Here a;-*- denotes the plane orthogonal to a; e , and we have used 
the fact that double-covers 63^2 homogeneously. 

Note also that by virtue of Example 3.6 and (3.9.1), we have 

(3.11.1) 7^;/3 = 7^;7^20 , 

where is the 1-form on given by (f){x) — 7r~^/^ e~l^'^v . 

Write r — \x\ , u — x/r , parametrize the unit circle perpendicular to 
u with unit speed as u^{9) , and put "w-polar" coordinates on by 
setting 

Then = r^\u — {u ■ u)u\'^ — r^(l — A^) , and hence 

n*f3{x) 

•1 r2Tr 



J^'e^"'" (2ttv- j\vLj^,0)u;x,e) dO^ dX , 



47r 

-r 



27r 



thanks to the symmetry of the integrand around A = . 

Now split V as v — v^^ + , with v^^ parallel to u , and • u — . 
A routine exercise then finds the parenthesized part of the integrand 
above to be 

TT (l + A") + 27r(l - A") . 

Put this into the last integral for 7^2/^ above, and make the substitution 
s :— r X to get 

7^2/?(x) = ll+{r)v^ + I_{r)v\\ , 
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where 

X±(r) = r-^e"'-' T e^' (r' ± s') ds . 
Jo 

Both X+ and X_ are smooth and decrease to zero as r — > oo . Specif- 
ically, it is not hard to show, using L'Hospital's rule, that 

lim r^X+(r) = 1 

r— »oo 

lim r^X_(r) = 1/2 , 
Putting these facts together with (3.11.1), we now see that 

ni(5 = 7^;7^2(/) ~7r( — + asr^oo 

In particular, we get decay, but not rapid decay, despite the rapid decay 
of the input field (5 . This is typical, as will follow from the Convolution 
formula (Prop. 4.4), and in this regard, the dual transform TZ\ behaves 
quite differently than the forward transform TZk ■ 

4. Convolution formula 

We now want to show that the composition 7^^ o 7^^ , applied to any 
locally intcgrablc and suitably decaying p-form on R" , convolves that 
p-ioim. with |a;|~^ 11 , up to some dimensional constant. As noted in our 
introduction, this fact generahzes the classical result of Fuglede [Fu]. 

Recall (§2.2) that 11 : R" — > Hom(/\^ R*^) is the operator field we get 
on R" by mapping each non-zero x e R" to the orthogonal projection 
IIj. : R" f\^Hx , where is the hyperplane perpendicular to x . 

The grassmannian Gn,k and the unit sphere S*^^^ C R" are both 
homogeneous spaces of 0(n) . To rewrite TZloTZ^ as a convolution, 
we will first pull integrals over the grassmannian back to integrals over 
0{n) , and then push them down to S"~^ . The following lemmas let 
us to do so with precision. 

4.1. Two averaging lemmas. The orthogonal group 0(n) acts on 
R" (and thus also on G„,fc) by left-multiplication. For any a G S"^^ , 
denote the a-stabilizing subgroup of 0(n) by Ka ~ 0(n — 1) . Then 
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0(n) is foliated by the mutually isometric left cosets of , and if 
7 e 0{n) , with b :— ja , we define the coset 

(4.1.1) Kab -.^iKa^ Kb-i = {« e 0{n) : Ka^b} . 

Note that 0(n) has a left-invariant metric whose Hausdorff measure 
assigns total mass 1 to both 0(n) and Ka — hence also to each coset 
Kab — in the appropriate dimensions. For integration, we always use 
these measures. 

We have adapted the first formula below from Helgason [H, I (15)]. 
We need the second one to handle differential forms, as opposed to the 
scalar-valued functions treated there. 

Lemma 4.2. Suppose maps Gn,k continuously into a vectorspace 
V . Then for any fixed Q e Gn,k , we have 



[ 4>{P) dP= [ 4>{iQ) 

JG„..k Join) 



in,k JO{n) 

Further, if a & S"""^ and %l) : 0{n) ^ V is continuous, then 
I ^i^) = I ( [ ^(«^) d^^) db . 

Jo{n) |S" 'I Js"-i KJk,, J 

Here db and dn denote the Hausdorff measures on S"~^ and Ka 
respectively. 



n,k 



Proof. Choose an arbitary /c-plane P G Gn.k , and map 0(n) — > G 
by defining ep(7) := 7P. Then ep{g'~^) = g ep{'y) for all g,jE 0(n) , 
and hence cp pushes the invariant probability measure on 0(n) down 
to the one on Gn,k ■ This lets us integrate over Gn,k relative to the 
latter measure by integrating o ep relative to the former. The first 
formula asserts nothing more than that. 

To get the second formula, we similarly map 0(n) S"~^ via 

60(7) 70- Then for any b e S"""-*^ , we have e~^{b) — Kab, and 
the coarea formula [Fe, 3.2.22] yields 



/ 

Joi 



^(7) Jca d-f = i ip{K) dn] d-f , 

'O(n) JS"-i \JK,i J 

with Jca denoting the Jacobian of . But this Jacobian must be 
constant, because commutes with the transitive left 0(n) action 
on itself. Set ip = 1 to see that Jca = \S^~^\ . □ 
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Notation: We write Ch(n, m) as a horizontal alternative to the usual 
vertical symbol for binomial coefficients: 

Ch(n,m) := (] = ^' rj . 

\mj ml[n — my. 

Lemma 4.3. Suppose P e Gn,k , x & P' and is the hyperplane 
perpendicular to x (of. ^2.2). Then for any a e Ap-^^" ^^^^ P , 
we have 

I Ch(fc,p) I 



/ 

J K.. 



Proof. Let A denote the "averaging" operator a ^ j^^ a\^p dn on 
the left-hand side of the formula above. Restriction to a subspace 
(followed by inclusion in /\pR"; see §2.2 (1) above) is a symmetric 
endomorphism of /\p R" , so A is symmetric too. 

Now observe that the action of on respects the splitting 

(2.2.3), namely 

It follows that A commutes with this action. But is the full 
orthogonal group of . So it acts transitively on both p- and (p — 1)- 
dimcnsional subspaces of , hence irreducibly on both l\pHx and 
X A /\p_iHx . This makes each summand in the splitting above an 
eigenspace of ^4 , so that in the notation of §2.2, we have 

(4.3.1) A^XU^ + ii^^ 

for some A, e R . But the definition of A makes clear that /jl — , 
which, modulo the particular value of A , proves our lemma. 

To evaluate A , observe that for any subspace V C R" , the projection 
Ap R-" — Ap^ trace given by 

tr(.|,)=dimA,V=(d';'^) 

on Apf^"- Taking the trace on both sides of (4.3.1) then gives 

tr(A) = A 

V P 

At the same time the definition of A gives 
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These two facts determine A and prove the lemma. □ 

We can now present the main result of this section: 

Proposition 4.4 (Convolution formula). Given a locally integrable p- 
form a on with < p < k < n and ||q;||^^) < oo , the convolution 
(r^'^n) T*r a exists, and we have 



S*^-^|Ch(A;',p) 



^>^'°= is-^-ich(;-'i:rt ''- 

Proof. Our assumptions ensure that for each fixed x G R" , the p-form 

y ' — ' \y\~^'^ya{y - x) 

is locally integrable and 0{\y\~"'~^) near infinity. This guarantees 
existence of the convolution. 

To see that the convolution encodes TZlTZk as claimed, combine the 
definitions of TZk and TZ^ with the fact that for any x G R" and 
P G G„ , we have Px + P' = x + P' as affine subspaces. It follows 
easily that 

(ni o 7^fe) a(x) = a(x + r]) |^ drj dP . 

JG„,k JP' 

Now fix an arbitrary Q G Gn,k and apply the first formula in Lemma 
4.2, then Pubini. Since ^Q' — {^r]: r] G Q'} for any 7 G 0(n) , this 
gives 

(4.4.1) (7^*o7^,)a(a:) = / / a{x + jrj)\ dj drj . 

Jq' J0{n) 

Note that when a has degree p = , it takes mere scalar values, and 
the restriction to 7(5 has no effect. In that case, the inner integral 
reduces to a simple x-centered spherical average of a , and the convo- 
lution formula quickly follows (cf. [H, Ch.I (34)]). Because we assume 
p > , however, the restriction complicates our task and we must now 
sort that out. 

Start by using the second identity from Lemma 4.2. If wc define 
fj = r]/\r]\ , and recall from (4.1.1) that Kfjgi = Kg/j when 7 G 0(n) 
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satisfies ^fj — q' , it rewrites (4.4.1) as 

TZlTZka{x) = I I 1/ a{x + \r)\q')\ dK] dq' drj 

l» \ Jq' Js"-i \Jk^,, ^ J 

= L (4 ^ '''' '''' 

where we have fixed some 7 G -fCf^g' . Since 7(5 is a /c-dimensional 
subspace of T^/S""^ , Lemma 4.3 now applies to give 

ni TZk a{x) 

= 777— rr / / if ctix + \r]\Q')\ ^ rffi; ) dg' dr; 

Ch(fc-,p) /■ /■ . ^1 I Ml w , 

ls^'-^lch(fc^p) / I fe'-i,,, 

- |S-|Ch(n-l,p) io is^-.^^^ + ^^^W^ 
switching from rectangular to polar coordinates on Q' , with r := |?7| . 
But by writing r*^'~^ = 7-"^! 7-^^^ |/ := —rq', and noting that = 
Hq' , we can make the reverse move on , to deduce 

^^^^"^"^^ |S-^|Ch(n-l,p) y^„"("-^)k^ 
Now simply recall from the end of §2.2 that restriction to Hy induces 
the projection we call Uy , so that 

a{x-y) r~'' = \y\~'' Uy a{x - y) . 

Hy 

Putting this into the equation above, we obtain our convolution for- 
mula. 

□ 

Remark 4.5. Though we don't include the proof here, the convolution 
above must decay near infinity. To be precise, one can show that as 
r := Ixl — >• 00 , 



|(r-'=r)*/(x)| 



O {r-') , ll/IUi < 00 
0{r-nnr) , = 00 
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for any focally integrable, vector- valued mapping / on R" , assuming 
that k <n, < oo , and T is a bounded operator field. 

Remark 4.6. As noted in Remark 3.3, our transform reduces to the 
classical Radon A;'-plane transform when p = . The reader will easily 
check that when we apply our convolution formula to any continuous 
0-form / : R" — > R with ||/||(„) < oo , it too reduces to a known 
result (cf. [H, p. 29 (55)]), namely 

Remcirk 4.7. Roughly speaking, TVjJZkOc will be at least as smooth 
as a itself. For when a has continuous derivatives D'^a for all multi- 
indices I7I <q, with \D'^a\ — 0{r~'^) , one easily checks that 

{{r-^ n) ^ a) = {r-^ Yi) ★ D^a , 

so that by our convolution formula, TZ^TZkCt is difi^erentiable through 
order at least q too. 

5. Kernels 

To fully exploit the convolution formula, we also need to know the 
Fourier transform of its kernel r~'' U . We assume k < n, so the 
singularity at r = is integrable, and r~'^ IT represents a tempered 
distribution. It therefore has a Fourier transform. Explicitly, 

Theorem 5.1. When < p < k < n , we have 

T (r-f'Tl) = 1 1^ {k-p)U+{n-p)^ 
^ ^ k jS''^'"^! r'^' 

When p < k , this makes T (r~*^n) pointwise invertible, with 

V J jS'^"^! \k — p n — pJ 

To prove this, we require two main ingredients. The first is Lemma 5.2 
below, a beautiful formula that Stein bases on an identity he credits to 
Hecke [St, p. 73]. The second. Lemma 5.5, prepares us to exploit Stein's 
lemma by expanding H as a sum involving spherical harmonics. 
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Lemma 5.2. Suppose : — > C is a homogeneous harmonic poly- 
nomial of degree d , and < k < n . Then r~'^~^hii and its Fourier 
transform are both tempered distributions, with 

Proof Stein [St, p. 73]. □ 

Remcirk 5.3. When d — and h = 1 , Lemma 5.2 makes the well- 
known classical assertion 

|S^'-V(r-^) = |S^-V~''- 

When p = , Theorem 5.1 above reduces to precisely the same asser- 
tion, because the projections 11 and $ reduce to the identity and zero 
operators respectively in that case, as noted in §2. 

Our next lemma handles the non-triviality of U and $ when p > . 
Recall that they are the operator fields on R" which, at each non- 
zero X G R" , project A^R" to /\^Hx and x/\/\^H,j. respectively (cf. 
2.2 A). The wedge operation and its adjoint (2.2.2) yield an elegant 
and useful pair of formulae for these operations: 

Observation 5.4. When 7^ x e R" , and a e A*^^"' > have 

\x\'^Iix{ci) = {x A a) y X 
\x\'^ $a;(Q;) = X A{a\J x) . 

Proof. Extend ei := x/\x\ to an orthonormal basis {ei,...,e„} for 
R". It then suffices to check both formulae on the corresponding mul- 
tivector basis ex defined by (2.2.1), an easy exercise. □ 

We call an operator field H : — > Hom(/\^R") harmonic if 
for all constant elements a,/? e A*!^"- 

Lemma 5.5. Let < p < n . Then there exists a homogeneously 
quadratic and harmonic operator field H such that, as operators on 
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p-forms, we have 

n = (p' /n)I + r-^H 
$ = {p/n)I-r^^H 

Here I is the identity on /\pR^ and p' = n — p . 

Proof. Since 11 + $ = 7 , it suffices to prove the first identity above, or 
equivalently, that 

nr^U^p'r'^ I + nH 

with H harmonic. 

Observation 5.4 says, for any a e Ap^^" > that 

nr^ n(a) = n (,t A a) V x , 

showing that nr^H is polynomially quadratic. Since the ring of ho- 
mogeneous quadratic polynomials on R" splits as r^C + , where 
Ti^ denotes the harmonic subspace, we can write 

(5.5.1) nr'^U = r'^A + nH , 

with A constant, and H e . Apply the laplacian to this expansion 
and divide by n to get 

(5.5.2) A(r^n)=2A. 

On the other hand. Observation 5.4 gives, for any basis p- vector of the 
form ex := ex, A Cx^ A ■ ■ ■ A ex^ (cf. 2.2.1) 

r"^ ^x{ex) = ^ XiXj (ci A ex) V ej . 



Since A{xiXj) = 2 6ij , and 



{ei A ex) V ej = 



ex, i ^ X 
0, i e X , 



we see that 

A (r^ n(eA)) = 2 ^(e^ A e^) V = 2 ^ = 2p' ca , 

Since the ex form a basis for /\pR,^ , this means A(r^ 11) = 2p' I , and 
comparison with (5.5.2) then yields A = p' I , turning (5.5.1) into the 
identity we sought. □ 
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Proof of Thm. 5.1. To compute the Fourier transform of r~'^n , first 
expand using the lemma just proven, to get 



Since / and H are harmonic, and homogeneous of degree and 2 
respectively. Lemma 5.2 computes the Fourier transform of this expan- 
sion as 



p' \S^-^\ I |S^+i| H 



(5.5.3) T{r n) - - ^ |gfc'+i| ^:Ff2 • 

To simplify this, recall that in terms of the classical Gamma function, 
one has 

lcfe+i| _ ^ ^ 



r(| + i) 

for any dimension k . But T{s + 1) = •sr(s), implying the recursive 
identity 

|cfc+i| _ ^ |S^~^| 
' ' A; ' ' ' 



which reduces (5.5.3) to 

|gfc'-i| ^k' \ k r'^ 



(5.5.4) ^(r-'n) = O-lf?:/-^! 



Now use the first identity in Lemma 5.5 to rewrite H/r'^ in terms of 
n and / , to get 



k' H 
k 


n 






l( 
n \ 


■4) 









n k r'^ n k \ n 



n' h' 

- |(n + *)-^n 

(A; - p) n + (n - p) $ 
k ' 

since p' — k' = k — p . Inserting this into (5.5.4), we obtain the desired 
formula for J^{r~'^lV) . 
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The subsequent formula for its inverse now follows immediately when 
p < k , because U and $ project onto orthogonal complements, so 



6. Inversion, Superposition, and Reconstruction from 



Finally, we will combine the Convolution formula (Prop. 4.4) with our 
computation of ^(r~'^n) (Thm. 5.1) to invert TZ^ exphcitly on the 
space of Schwartz p-forms. 

When p — , our formula reduces to the well-known inversion of the 
classical /c'-plane transform. The latter result apphes well beyond the 

Schwartz class, to all continuous functions [H, Thm. 6.2] (or even func- 
tions in L«,^, with 1 < g < n/k' [R, Thm. 4.1]), with 0(r-"-^) 
decay. 

We expect that, likewise, our inversion formula holds for p-iorms with 
0{r~"'~'^) decay, leading to both a superposition formula, and an algo- 
rithm for reconstructing currents from their projections. We will state 
the results we anticipate along those lines. We leave their proofs, how- 
ever, to colleagues more expert in the analytical techniques they seem 
to require. 

Theorem 6.1 (Inversion). Suppose k < n, and a is a Schwartz p- 
form on R" , with p < k . Then we may recover a from its transform 
TZka as follows: 



Here T and T ^ apply in the tempered distributional sense, and 



that n$ = $n = and R + $ = 7 . 



□ 



Projections 




'n,k,p 



k 



S"-i|Ch(n- l,j9) 
S'=-i|Ch(n- A;,p) 



Proof. Apply the Convolution formula (Prop. 4.4) to get 

nlTZka^ ci (r-'^R) ★ a , 



where 



S'='-i|Ch(A;',p) 



S"-i|Ch(n-l,p) ■ 
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As noted at the start of §5, H represents a tempered Hom(/\^R'*)- 
valued distribution on R". The multiplication rule for T therefore 
extends to the convolution above (see Observation 7.2 for a proof), 
and we can write its Fourier transform as 

Theorem 5.1 computes J^ir^^Ii) explicitly as an operator field, and 
shows that the right-hand side above, though apriori just a tempered 
distribution, is actually represented by a locally integrable p-form. We 
can then multiply both sides by the field pointwise inverse to that p- 
form (also given by Theorem 5.1) to deduce 

^(«) = ^-^T^' (-^-r^\ .T{ninua) 

Cl \k — p n — p J 

where 

|gfc'-i| 

(6.1.1) c2^C2{n,k) -.^k-^-j^ . 

The Fourier transform preserves Schwartz spaces, so T{a) , fike a 
itself, belongs to SpiW") . In particular, the right side of the formula 
above lies in <Sp(R") , and we can therefore recover a from TZ^a by 
applying the inverse Fourier transform: 

(6.1.2) a=''^T-'(r^' + • ^ (7^^ 7^,a)') . 

Cl \ \k-p n-pj ) 

Since C2/C1 = kCn^k,p, this completes the proof. □ 

When the codimension k' = n — k is even, the unwieldy formula above 
takes a simpler "differential" form, using the operators d and 5 instead 
of the Fourier transform: 

Corollary 6.2. When the hypotheses of Theorem 6.1 hold with codi- 
mension k' —: 2Z e 2Z, we can recover a from TZ^a via 

Proof. By the formulae in Proposition 2.9, JF intertwines d5 and 5d 
with multiplication by and r^II respectively. But $ and 11 are 
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projections, so = $ and 11' = 11 for any integer Z > , whereby 

When k' — 21 , these facts simphfy our inversion formula in an obvious 
way. Remark 4.7 then ensures that TZlTZkOi is suitably differentiable, 
and the corollary follows. □ 

Remark 6.3. The inversion and superposition results above all cor- 
rectly reproduce their classical scalar precedents when p = , since 11 
is the identity and $ = on /\q(R") . Indeed, one easily checks that 
these facts reduce the inversion formula of Theorem 6.1 to 

for any Schwartz function / . For functions, 6d = and d5 is the 
Laplace operator A . So when k' = 21 is even and p = , Corollary 
6.2 reduces to 

in agreement with [H, Thm. 6.2]. 

Indeed, when p = the inversion formula extends beyond iSp(R") to 
(for instance) the class of all continuous functions a with 0(r~"~^) 
decay near infinity, as mentioned above, and we expect the analogous 
extension to hold for p > : 

Conjecture 6.4. The inversion formulae of Theorem 6.1 and Corol- 
lary 6.2 hold for all continuous p-forms a with 0{r~^^^) decay near 
infinity for some £ > . 

Now recall again that when a e »Sp(R") , so is J^a . This fact makes 
it easy to show that 

(6.4.1) ( r'' ( + . Ta\ 

\ \k — p n — p J J 
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is smooth, with an 0{r decay estimate. In particular, is inte- 
grable on any /c'-plane, so that 



(with Cn,k,p as in Theorem 6.1) is weU-defined and smooth on Tn,k ■ A 
shghtly better 0(r~"~^) decay bound for would combine with our 
Conjecture above and our inversion formula to let us write 



is clearly invertible, and hitting both sides of (6.4.2) with its inverse, 
given the definition of in (6.4.1), would then give a = 7^1(3 , making 
a a superposition of /c-planar p-forms. So if Conjecture 6.4 holds, so 
does 

Conjecture 6.5 (Superposition). Given < p < k < n , any a e 
<Sp(R") can be constructed as a superposition of k-planar p-forms by a 
formula of the type 



where (f) is given by (6.4-1) above. 

As we have noted, this has been proven when p = ([H, Thm. 6.2]). In 
view of Remark 3.10, it then analyzes the scalar- valued Schwartz func- 
tion q; as a superposition of "/c'-plane waves," i.e., functions invariant 
with respect to the translations generated by k' independent vectors. 
That fact accounts for the classical use of Radon transform methods 
to solve certain linear partial differential equations by the "method of 
spherical means." (cf. [H, §7a] and [R, (1.1)]). 

Conjecture 6.5 would similarly let us analyze Schwarz p-forms as super- 
positions of /c-planar forms. This in turn would allow us to reconstruct 
p- currents from their projections. 

To explain this, we introduce the linear space (R-"^) comprising all 
smooth p-forms on R" . One topologizes (R"^) in a standard way, so 
that a sequence in (R") converges to zero if and only if its members. 



P '■— Cn,k,pT^k4> 



(6.4.2) 




But the operator 




« = Cn,k,pT^lT^k(t> , 
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together with all their derivatives, converge uniformly to zero on any 
fixed compact subset of R" (cf. [H, Chap. V]). 

Definition 6.6 (Currents). We call any continuous linear functional 
on (R") a compactly supported p-current on R" . 

One easily shows that any such current actually does have compact 
support, in the sense that for some i? > , it vanishes on all smooth 
p-forms supported outside the ball Br{Q) . Since we can find a Schwarz 
p-form that agrees, on Pij(O) with any given smooth p-form, it follows 
immediately that 

Observation 6.7. A compactly supported p-current is determined by 
its values on 5p(R") C £p (R") . 

Finally, note that wc can map any such current T into a A;-plane 
P G G„,fc , to produce a compactly supported A';-currcnt called the pro- 
jection of T into P , and denoted by P*T, via the simple rule 



Here, as in our introduction, P*(f) denotes the puUback of (f) from P 
to R" via the orthogonal projection R" — > P . 

Conjecture 6.8. When p < k < n , each compactly supported p- 
current T on R" is uniquely determined by its projections into all 
k-planes P e Gn,k ■ Explicitly, for any a e <Sp(R"') , we have 



J Gn,k 

where (5 is the mapping constructed from a in Theorem 6. 5. 

Proof (modulo Conjecture 6.4)- By Observation 6.7, it suffices to show 
that whenever a e (Sp(R") , we can determine T[a] from the projec- 
tions P*T for every P G G„ ^ . But if Conjecture 6.4 holds, then so 
does Conjecture 6.5, and the latter constructs a p-form from a for 
which the mapping /3 := TZk4> satisfies 



P,r[0] r[P>] for all G £^ (P) . 
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for every x e R" . Since (3 lies in the image of TZk , Remark 3.4 makes 
the restriction |p superfluous in the formula above, which we can thus 
rewrite as 



To clarify the main logic of our exposition, we have deferred two tech- 
nical facts about convolution to this appendix. 

We start by verifying that the notion of convolution extends to the 
setting required by our inversion theorem. 

Observation 7.1. When F G 5(R", V") , the classical definition of 
T-kF (see Item (3) in %2.7) extends in a natural way to the case where 
T is a tempered {V^W)- distribution. 

Proof. To do this, we identify vectors with operators in a somewhat 
unusual way. 

Let V* denote the image ol v &V under the canonical duahty V ~ 
V* . We then let v^:W ^ Hom(y, W) signify the operator we obtain 
by setting 

{v^{w)) {x) := {w ® v*){x) := v*{x) w eW 

for all a; G and w & W . We can then identify any y-field F on 
R" with the field F^ given by 

X ^ F^{x) := F{x)^ . 




It would then follow that 




proving the result. 



□ 



7. Appendix: Technical lemmas 
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A routine calculation shows that whenever T e »S(R"', V, W) , F e 
5(R", V) and G e 5(R", W^) , we have 

(7.1.1) {T . F)[G] = T[G ® F*] . 

Using this fact with a bit more calculation then reveals that 

(7.1.2) (r*F)[G]=r[F-*G] , 

where the superscript "— " indicates reflection through the origin, e.g., 

F-{x) := F{-x) . 

Since F^ i< G lies in 5(R", VT) (7.1.2) now deflnes T i< F for any 
T e cS* V, W) , as promised. □ 

We next show that the Fourier convolution rule (Item (3) of §2.7 again) 
extends to convolutions of the type just described: 

Observation 7.2. The Fourier convolution/product rule 

T{T-kF)^ T{T).T{F) 
holds when T e 5*(R", V, W) and F e 5(R", V) . 

Proof. In view of (7.1.1) and (2.7.2), we have to show that for any 
G e 5(R", W) , 

(7.2.1) (T ★ F) [T{G)] = T [r{G ® T{F)*)\ . 

We start by using (7.1.2) to write 

{T^F)[T{G)]^T[F-^T{G)\ . 

Since F and G are both Schwartz, properties (1), (2) and (especially) 
the convolution/product rule (3) listed above, together show that 

F-^^(G) = t-'{t{f-).:f\g)) 
= T{G®T{Fy)) . 

The desired fact (7.2.1) follows immediately. □ 
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